NACA RM L56 f J15 



Copy 339 

BM L56J15 



RESEARCH MEMORANDUM 


APPLICATIONS OF POWER SPECTRAL ANALYSIS METHODS TO 

MANEUVER LOADS OBTAINED ON JET FIGHTER 

AIRPLANES DURING SERVICE OPERATIONS 

By John P. Mayer and Harold A. Hamer 

Langley Aeronautical Laboratory 
Langley Field, Va. 


CLASSIFIED DOCUMENT 

Ttds material contains Information affecting the National Defense of the United States within the 
of the espionage laws, Title 18, U.3.C., Secs. 793 and 794, the transmission or revelation of which in any 
manner to an unauthorised person la prohibited by law. 

NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

January 15, 1957 



TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 


NACA EM L56J15 

NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


limit 

014422S 


RESEARCH MEMORANDUM 


APPLICATIONS OF POWER SPECTRAL ANALYSIS METHODS TO 
MANEUVER LOADS OBTAINED ON JET FIGHTER 
AIRPLANES DURING SERVICE OPERATIONS 
By John P. Mayer and Harold A. Hamer 

SUMMARY 


Power spectral densities of normal load factor have "been obtained 
for two service operational training flights of a Republic F-84G air- 
plane and three service operational training flights of a North American 
F-86A airplane in order to indicate the load-factor frequency content 
and possible uses of power spectral methods in analyzing maneuver load 
data • 


It was determined that the maneuvering load-factor time histories 
appeared to be described by a truncated normal distribution. 

The power spectral densities obtained were relatively level at fre- 
quencies below 0.03 cycle per second and varied inversely with approxi- 
mately the cube of the frequency at the higher frequencies . In general, 
the frequency content was very low above 0.2 cycle per second. 

The load-factor peak distributions were esti ma ted fairly well from, 
the spectrum analysis. In addition, peak load data obtained during 
service operations of fighter-type airplanes with flight t ime totaling 
about 2^,000 hours were examined and appeared to agree reasonably well 
with the type of equations obtained from spectrum peak-load distributions. 


INTRODUCTION 


With a view toward providing Information for the appraisal of the 
design requirements of future airplanes, the National Advisory Committee 
for Aeronautics with the cooperation of the U. S. Air Force and the 
Bureau of Aeronautics, Department of the Navy, has been conducting a 
statistical study of the pilot input, the airplane motions, and the 
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loads imposed on a number of fighter-type airplanes during military 
service operations. The data obtained thus far are - summarized in 
reference 1. 

The statistical treatment of the maneuver loads has been limited 
to the peak values of the loads and the frequency of' occurrence. In 
the analysis of other loads associated with such random factors as gusts, 
buffeting, or runway roughness (for example, refs. 2 to 5)> the applica- 
tion of power spectral analysis methods has proved valuable, especially 
since a number of useful theoretical relationships such as the input- 
output relations and the frequency of occurrence of peak loads have been 
established for stationary random processes. Such methods of analysis, 
however, have not been applied to maneuver loads . A single maneuver is 
clearly not a random process since It is associated with pilot intent. 
Nevertheless, maneuvering of a general nature (other than specific 
missions) over a long period of time or a collection of a great many 
maneuvers may acquire same of the characteristics o£_a random process. 
Therefore, it was believed that some of the relationships from the theo- 
ries of spectral analysis might be used as a guide in analyzing the 
maneuver load data. 

With this in mind, the normal load factors obtained in several 
flights of two jet fighter airplanes flown In regular squadron operations 
were analyzed by using power spectral methods. The data analyzed were 
obtained during the investigations reported in reference 1 and they are 
a portion of the results summarized in reference 1. 

It Is the purpose of this paper to present the power spectra obtained 
in order to show the frequency content of normal load factor in service 
operations and to determine whether other information obtainable from 
spectral analysis methods (such as the frequency of occurrence of peak 
loads) could be obtained from the maneuver load spectrum. 

Power spectral densities at frequencies up to 1 cycle per second 
are presented for two operational flights of a Republic F-84G airplane 
and for three operational flights of a North American F-86A airplane. 

The peak loads derived from these spectra are compared with the actual 
peak loads. In addition, peak-load data from other service airplanes 
are examined and compared with the results of spectrum analysis. 


SYMBOLS 


a constant used in spectrum equation 

f frequency, cps 
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-B 


"break" frequency, cps 

frequency of occurrence of normal load factor 

average number of normal load factor crossings per second 
of the zero axis with positive slope. 


*o 2 = 


Jp 

poo 

L * (f) ^ 


df 


(ref. k) 


df 


total number of positive normal- load-factor peaks per second. 
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(ref. 4) 


m 

N 

n 

n(t) 

n L 

An 

An L 

An 

An" 

Anp 

P(y) 


total number of observations 

total number of positive peaks exceeding a given value 

normal load factor, g units 

time history of normal load factor 

service limit normal load factor 

incremental normal load factor (n - 1) 

incremental service limit normal load factor ^n^ - 1^ 

mean incremental normal load factor 

filtered value of incremental normal load factor; at time t, 

= °-9^a fc _At “ S£ * L t + °-9^t+At 

value of An at point of truncation 

probability that a given value y will be exceeded 
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P N (y) 


T 

% 

t 

At 

a(0 


8(f-0) 


a 




®(fW" 

X 2 



probability that a given value y w ill be exceeded, based on 
a normal probability distribution 

probability of exceeding a given load-factor peak for load 
factors greater than 1 

total flight time ’ 

maneuvering flight time 

time, sec 

time reading interval, sec 

proportion of area of normal curve represented by data in a 
truncated distribution 

Dirac delta function at f = 0 — - — 

mean of parent distribution from which truncated distribution 
is derived 


root -me an- square value of 
power spectral density of 


0 z(An - An) 2 r 00 . 

^1, a 2 = = / 

m lo 

incremental normal load factor (An), 



cps 


> 



(ref. 2) 


power spectral density of An" 
statistical reliability parameter 
critical statistical reliability parameter 


AIRPLANES 


The airplanes for which data were obtained were service models of 
the Republic F-84G and the North American F-86A. Both are low-wing 
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Jet-propelled fighter-type airplanes, the F-86A having a swept wing and 
empennage . The F-86A was equipped with an adjustable stabilizer and a 
hydraulically boosted elevator. 

The power spectra presented herein were obtained from flights made 
with external fuel tanks for the F-84G and without external fuel tanks 
for the F-86A. Except for the addition of sideslip and angle -of -attack 
booms, neither the external appearance nor the weight and balance of the 
airplanes was altered by the addition of the NACA instrumentation. 

Three -view drawings of the airplanes are presented in figure 1. 


INSTRUMENTATION AND TESTS 


The data were measured by standard NACA photographically recording 
instruments. In order to relieve the pilot of any recording-instrument 
switching procedure and thus to assist in obtaining normal operation, 
a pressure switch was employed to operate the recording instruments auto- 
matically at take-off. 

Normal load factors (normal to the longitudinal reference) were 
measured by an NACA air-damped recording accelerometer. The measuring 
element was damped to about O .65 of critical damping and the natural 
frequency was about 20 cycles per second. The accelerometer was located 
near the airplane center of gravity in such a way that the effects of 
angular velocities and angular accelerations were negligible. 

A standard two-cell pressure recorder connected to the airplane 
service system was used to measure the pressure altitude and indicated 
airspeed. The service systems were of the usual total-pressure -tube and 
flush static-pressure-orifice type. No corrections were made to the 
measurements for position error associated with the location of the 
static orifices. 

The flights were performed by service pilots during regular squadron 
operational training. Data were recorded continuously throughout a 
flight and were recorded only during those flights in which the mission 
was scheduled to include a large number of maneuvers. Although not 
requested, most of the maneuvers were performed in relatively smooth air. 

No attempt was made to specify the type or severity of maneuvers. Although 
the pilots were aware of the instrumentation. It was stressed that this 
wan not to restrict their normal handling of the airplane since they would 
not be personally identified with the test results. 

The flights for which power spectra were obtained were selected since 
they contained most of the tactical maneuvers that are within the capa- 
bilities of the individual airplanes. A summary of the operations and 
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flight conditions for these flights is presented in the following table 
and includes the actual amounts of time from which the spectra were 
obtained. Distributions of the percentage of time (used for spectra} 
spent in various altitude and airspeed ranges for two combined flights 
of the F-84G and three combined flights of the F-86A are shown in fig- 
ures 2 and 3, respectively. 


Airplane 

Flight 

Pilot 

Ml 68 Ion 

Average 

weight, 

lb 

Average 

pressure 

altitude, 

ft 

Average 

indicated 

airspeed, 

knots 

Average 

Mach 

number 

Time for which 
spectrum was 
obtained, sec 

Republic 

F-84G 

1 

A 

Acrobatics and 
dive-bombing 

16,000 

11,000 

290 " 

0.53 

3,900 

2 

B 

Acrobatics, 
dive-bombing, 
and ground 
gunnery 

2,777 


B 

C 

Acrobatics 





1,876 

North American 
F-86A 

2 

C 

Aerobatics 

13,060 

14,150 

293 

.58 

1,823 


3 

O 

Dogfighting 





1,032 


METHODS OF ANALYSIS 


Power Spectral Densities 

In the present paper use will be made of some of the concepts from 
the theories and methods of power spectral .analysis. The methods used 
in the present paper are described in references 2 and 3. 

The methods described in reference 2 were used in determining the 
power-spectral densities of the normal -load-factor data. The records 
were read at 0.5-second intervals for each flight and kO estimates of 
the power at frequencies up to 1 cycle per second were obtained by using 
digital computing methods. In order to improve the accuracy of the com- 
puted spectra, the spectra were first obtained of an altered (filtered) 
time history of load factor given by the equation 


( 1 ) 


An/ = 0.9 dnt_At - 2 An^. + 0.9 ^t +At 
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This is essentially a high-pass filter which retains some information 
at zero frequency. This particular expression was used since it was 
found to make the spectrum of in" relatively constant over the fre- 
quency range. The spectrum of in" was then converted to the spectrum 
of in hy the equation 


*(f) 


0(f) 


in" 


in , . 2 

(l.8 cos 2itf it - 2) 


( 2 ) 


In addition to the spectra computed from the records read at 
0.5-second intervals, spectra were also computed from the data taken 
at 2-second intervals and at 10-second intervals in order to provide 
more detail on estimates of power at the lower frequencies. The spec- 
tra computed for 2-second intervals were computed for each flight at 
frequencies up to l /4 cycle per second. The spectra computed for 
10-second intervals were computed for the combined individual flights 
for each airplane at frequencies up to l /20 cycle per second. The data 
at 10-second intervals were not filtered by using equation (l) since 
the spectrum was relatively level. 

Ordinarily spectra are obtained of quantities which fluctuate about 
some mean value. Maneuver load time ‘histories, on the other hand, fluc- 
tuate about 1 g but generally only in the positive direction. 

The power spectral densities of the incremental maneuvering load 
with a nonzero mean (An) would include a Dirac delta function at the 
origin with an area ( 55 ) 

= 2^8(f-0) + (5) 


The power spectra presented in this paper have the mean taken 
out and are equal to the last term ®(f )^ n in equation (5)'. 


Load-Factor Counts 

Threshold counts . - The probability curves for the normal-load- 
factor data were determined by a process which minimized the level- 
flight portions of the record. In this method counts were made at a 
given threshold every time the record crossed the given threshold and 
exceeded one -half of the increment to the next threshold. In this paper 
threshold values of 0 > 5 g were ■used. For example, counts were made at 
2g if the load-factor record crossed 2g while increasing and then 
increased beyond 2.25g or if the record crossed 2g while decreasing 
and decreased beyond 1 . 75 g> This method Is Illustrated in figure 4 
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(method 3)* In using this method all small fluctuations in the load 
factor less than +0.25g are eliminated. 

Peak counts . - Two methods of obtaining major peak distributions 
were used in this paper and are illustrated in figure k. In the first 
method (method l) peaks were counted with a threshold of 0.25An or 
0.25g whichever was greater; that is, the load factor must decrease _ 
more than one -quarter of the value of the peak -load-factor increment 
An before and after the peak. For values of Ai below 1 (2g), how- 
ever, a threshold of 0.25g was used in which the load factor had to 
decrease 0.25g before and after the peak. The second method (method 2) 
was method B of reference 1. In both methods peaks of less than 0.25g 
were not counted. Peaks obtained by using method 2 -were counted manu- 
ally. Those obtained by using method 1 were counted automatically from 
the time histories with a digital computer. 


RESULTS AND DISCUSSION 


Spectrum for F-8LG Airplane _ 

In figure 5 the power spectrum of normal load factor is shown for 
flight 1 for the F-84G airplane. Results -are shown for the four quar- 
ters of the flight, of which each part consisted of about 900 seconds 
of flight. The root -mean-square value cx and the mean Zn for each 
pant are shown in the figure. 

The power spectral densities for each of the four parts appear to 
be very similar and decrease approximately with f~ 5-. 

The total spectrum for flight 1 of the F-84G airplane is shown in 
figure 6. Also included is the spectrum of the data taken at 2-second 
intervals to provide more detail at the lower frequencies. In addition, 
the statistical 95 -percent confidence bands associated with the spectrum 
are shown. For simplicity, the confidence bands shown were obtained 
for a faired curve rather than for each individual point. The value of 
^(f )ah zero frequency shown has little statistical reliability and 
and is shown only as a matter of interest. 

The confidence bands shown in figure 6 are representative of those 
obtained for all the spectra presented in this paper. Therefore, con- 
fidence bands for the other spectra will not be given. 

The spectrum obtained with At = 2 seconds provides more detail 
at the lowest frequencies and generally agrees with the spectrum obtained 
with At = l/2 second except at the highest frequencies. At the highest 
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frequencies for the 2-second data (l.e. , l/h cps) the spectrum is not 
expected to he accurate since the power at higher frequencies is included 
in the spectrum obtained. Therefore the values for the 2-second data 
are not shown at the higher frequencies. 

Effect of level flight on spectrum . - The effect of the presence in 
the time history of long intervals at 1 g or level flight on the spec- 
trum is shown in figure 7* The results shown are for about 900 seconds 
from flight 2. The record was examined and approximately 550 seconds 
of flight near 1 g were deleted. The resulting record was then analyzed 
with one maneuver connected to the end of the previous one. The results 
are presented as the power spectral densities divided by the mean square 
<j2 0 f the respective parts. The value of 02 0 f the original was 0-598 
and the value of f or the all- maneuvering time history was 0.929* 

It may be seen from figure 7 that above 0.025 cycle per second the 
spectra are directly proportional to the mean square or 





L « 2 

T 

_ 02 


(4) 


where T^ refers to the maneuvering flight time and T to the total 
flight time. In addition. 


° 2 t _ % _ 


T 

-2 T 

° % 


TM 


(5) 


In general, the effect of level flight on the spectrum is to change the 
level of the spectrum. 

Spectra for flights 1 and 2 . - In figure 8 the spectra for the com- 
plete flights 1 and 2 are given. In each case the power has been divided 
by the mean square for the particular flight. The power spectral den- 
sities are seen to be quite similar for the two flights. At the higher 
frequencies (above 0.1 cps) the two curves have the same shape but are 
of different magnitudes . 

In figure 3, flights 1 and 2 have been combined and, in addition, 
the power spectral densities obtained for the combined flights 1 and 2 
at 10-second Intervals have been included. The 10-second data have been 
corrected for "folding" based on the spectra obtained from the l/2-second 
and 2-second data. Wien the 10-second data are included, a peak appears 
at about 0.025 cycle per second which corresponds to the average number 
of maneuvers per second during the two flights, which is roughly 0.024 
(86 maneuvers per hour). 
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Spectra for F-86A Airplane 

The power spectral densities for three operational flights of the 
F-86A airplane are shown in figure 10. The results are ratioed to the 
mean square values of the incremental load factor. The spectra for 
these three operational flights also appear to he similar. 

The power spectral densities for the combined flights 1, 2, and 3 
for the F-86 airplane are shown in figure 11. Included in this figure 
are the data obtained at 10-second intervals. As was the case for the 
F-84G airplane , there appears to be an indication of a peak in the spec 
trum at approximately 0.02 cycle per second. Again this value corre- 
sponds to the average number of maneuvers per second in flights 1, 2, 
and 3 of about 0.022 (79 maneuvers per hour). 


Comparison of Spectra for F- 8 kG and F- 86 A Airplanes 

In figure 12 are shown the normalized power spectral densities for 
three combined flights for the F- 86 A airplane and two combined flights 
of the F- 8 kG airplane. The same curves are shown on a linear scale in 
figure 13 . 

The power spectral densities for the two airplanes in squadron 
operations are very much alike; however, there appears to be some 
increase in power at the highest frequencies (l cps) for the F -86 air- 
plane as compared with the F-84G airplane. It is believed that this 
increase is caused by a lightly damped longitudinal oscillation which 
was present In the F -86 airplane. The effect of this is probably pres- 
ent at frequencies between 0A and 1 cycle per second. The power is 
very low above 0.2 cycle per second. (See fig. 13 .) 


Analytical Representation of Maneuver Load Spectrum 

The spectra of figure 12 appear to be of the type given by the 
formula for the simple spectrum: 


= 


1 + (f/f B ) a 


(6) 


where f B is the "break" frequency. In figure ih- the total spectra 
for the F- 8 LG and F- 86 A airplanes (shown in fig. 12) were averaged 
according* to the flight time for each airplane and are conpared with 
equation ( 6 ). In this case a was found to be 2.6 and f B was 0.051 
cycle per second. This equation appears to fit the data fairly well. 
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Also shown in the figure as a matter of interest is the approximate 
phugoid frequency range which is between 0.008 and 0.025 cycle per sec- 
ond. In addition, the average number of maneuvers per second was 0.022 
and 0.024 per second for these flights (79-86 per hour) . 


Probability Distributions 


‘Threshold distributions . - The frequency distribution of threshold 
counts of normal load factor for the combined three flights of the 
F- 86 A airplane is shown in figure 15 and table I. In examining the 
frequency distributions, it was noted that the distributions looked 
much like normal distributions if only the positive or only the nega- 
tive distributions were examined. For example, in figure 15 , normal- 
distributions curves are shown for both the positive and the negative 
distributions. These may be described as truncated normal distributions 
above and below An = 0 . (See ref. 6 for a discussion of truncated dis- 
tributions.) Since the distributions are essentially one sided, the 
measured value of the threshold counts at An = 0 must be adjusted. 
Therefore, in figure 15 it may be noted that the zero value for the 
positive side is given as two times the number of zero threshold values 
having a positive slope, and, in a similar manner, the zero value for 
the negative side is given as two times the number of zero threshold 
values having a negative slope. In table I the threshold counts- 
labeled 0 are the zero values having positive slope and those labeled (-)O 
are the zero values having negative slope. 


The probability of exceeding a given incremental normal load 
factor An-L for a truncated normal distribution is 


p (^i) 


1 

a (5)\/St 




dZ 


( 7 ) 


for An^ ?? 0 
where 


An^ = An - Auq 

z.^sni 

o 


and 



12 


NACA RM L 56 JI 5 


The normal -load -fact or data are truncated at Ac = 0 and the mean 
of the parent distribution is assumed to be at An = 0. Therefore, 

a( 0 = 0.50 

| = 0 

and 



Since it is sometimes convenient to examine data on normal proba- 
bility graph paper, the quantity a( £)P(An) of a truncated normal dis- 
tribution can be plotted and the results will fall in a straight line 
on normal probability paper. ' — 

The probability distributions of threshold counts of normal load 
factor for flights 1 and 2 of the F-84G airplane and flights 1 , 2 , and 3 
of the F- 86 A airplane are shown on a normal probability scale in fig- 
tires 16 and 17 . The data are shown for the individual flights as well 
as for the combined flights for each airplane. The value cc.(f;)P(An) 
is plotted against An. In each case, the zero value used was two times 
the number of zero thresholds having a positive slope. The threshold 
counts for each airplane are given in table I. The points shown in 
figures l 6 and 17 are plotted at l/4g less than the threshold values 
used (i.e., 1 plotted at 0.75^ 1-5 plotted at I. 25 , and so forth) since 
the thresholds were based on bands +l/4g about the threshold values. 

In figures 16 and 17 it may be seen that the truncated normal dis- 
tributions appear to fit the data. In addition to the plots on the 
normal probability scale the (Chi-square) statistical test was made 
for each flight and the combined flights. (See ref. _7, for example.) 

The results of these tests as shown In the figures are well within the. 
95-percent level of significance and indicate that the maneuvering loads 
of the operational flights presented in this paper can be considered as 
truncated normal distributions . 

Although similar analysis could be made of the negative maneuvering 
load factors, the data sample sizes used were not sufficient for any 
reliable analysis to be made. 

Peak-load distribution . - The peak -load -factor frequency distribu- 
tions obtained for the combined flights of. each of the airplanes are 
given in table I. Two peak distributions are given: the total peak 

distribution and the major peak distribution. The total peak distribu- 
tion includes all peaks obtainable from the time histories of normal 
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load factor. Since the load factors were calculated to O.OOlg, the 
threshold for the total distribution is O.OOlg, which is less than the 
reading accuracy of about O.Olg. A count was made of all the peaks by 
using a threshold value of O.Olg for the F-84G- airplane and it was found 
that the difference between the frequencies counted by using the two 
thresholds was almost entirely included in the 0 to Q.25g range. There- 
fore, the values beyond 0 . 25 g are probably not affected to any extent 
by the threshold used. In addition, since the records were read at 
0 . 5 -second intervals the minor peaks at high frequencies were also fil- 
tered out; however, these are not believed to be of any importance in 
these flights. 

In figures 18 and 19 the major peak cumulative frequency distribu- 
tions are shown plotted against (faa) 2 for the F-84G- and F- 86 A airplanes. 
Also plotted in figures 18 and 19 are the threshold frequency distribu- 
tions. From equation ( 8 ) it may be seen that the threshold frequencies 
should plot as a straight line when plotted against (An) 2 on semilog 
paper if the distribution is a truncated normal distribution. Also it 
may be seen from figures 18 and 19 that the major peak cumulative dis- 
tribution plots approximately as a straight line against (An) 2 above 
An = 1. 

Calculated peak -load distributions .- In reference 4- an expression 
for the probability of exceeding a given peak is given for a Gaussian 
random process. Although the maneuvering load data of this paper can- 
not be considered a Gaussian random process, the results of figures l 6 
and 17 indicate that the positive maneuver loads can be considered as a 
truncated or "half" normal distribution. Therefore, it was believed 
that perhaps the peak probability expression given in reference 4 might 
be used as a guide in calculating the peak distributions obtained in 
■the present tests. 


If the maneuvering loads were Gaussian in character the probability 
that a peak load will exceed a given value can be given, as in refer- 
ence 4, as 



where 


Kl = \/l - (fb/fp ) 2 


and 


K 2 = 


Kl 

f o/ f i 




>» «• 



Ik 
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and where the normal probability distribution Pjf(y) is 



where i = 1 or 2 . 

The total number of positive peaks exceeding a given value is 

N = f p %P(An) (11) 

For large values of An equation (9) approaches 

l/An\ 2 

P(Aq) « £o e ( 12 ) 

f p 

and for the asymptotic case the number of peak loads exceeding a given 
value approaches 


-ifef (15) 

N « f 0 %e 2 \ 0 > 

Since much of the flight time is spent in nonmaneuvering or near- 
level flight, the maneuvering flight time rather than the total flight 
ti m e is used in equations (ll) and (13). In order to determine the 
maneuvering flight time, the percentage of total flight time spent above 
a given load factor was determined. These times are shown in figure 20 
for the F-84G airplane and in figure 21 for. the F-86A airplane. The 
distribution of the time spent in various load-factor Intervals is given 
in table II. Since in the threshold counts and in the_peak counts the 
minimum load -factor variation was 0.25g, It was assumed that the maneu- 
vering flight time was the portion spent above An = 0 . 25g • It can be 
seen in figures 20 and 21 that about 52 percent of the.. total flight 
time could be considered maneuvering flight time for both the F-84G air- 
plane and the F-86A airplane in these particular flights. 

The peak -load data for the F-84G airplane are shown in figure 22 
and the peak-load data for the F-86A airplane are shown in figure 23. 

Two peak-load distributions are shown, - the total peak distribution 
and the major peak distribution. The major peaks shown^In figures 22 
and 23 were counted by both method 1 and method 2 of figure 4; however, 
only the results from method 1 are shown. The results from method 2 
were almost the same as those from method 1. 
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The peak-load distributions as given by equations (ll) and ( 13 ) 
are also shown in figures 22 and 23- The value of the root -mean-square 
load factor a used was that obtained from the threshold analysis of 
figures 16 and 17 . The value of ^oAp obtained from the integration 
of the spectrums was 0.184 for the F-84G airplane and 0.135 for the 
F-86A airplane. The value of fQ was O.O 76 for the F-84G airplane and 
0.081 for the F- 86 A airplane. 


It may be seen in figures 22 and 23 that both the total peak count 
and the major peak count are approximately estimated by equations (ll) 
and ( 13 ) , respectively. The agreement seems somewhat better for the 
F- 86 A airplane than for the F-84G airplane. 


Although the total peak count is of interest, it is the major peak 
distribution which is of most value. For example, in fatigue analysis, 
it is believed that the major peak distribution would more accurately 
describe the life history of loads than would the total peak distributions 


Major peak -load distribution in service operations .- In order to 
compare the results from the spectrum analysis with more complete flight 
data from regular service operations the major load-factor peaks obtained 
in the service operations of references 1 and 8 to 11 are shown in fig- 
ures 24 to 26 . The airplanes for which information Ib presented in 
figures 24 to 26 are 


Republic F-84E and G 

North American F- 86 A, E, and F 

North American T-28A 

North American B-45A 

North American AJ-1 

Lockheed F-80A, B, and C 

Lockheed F-94B 


Lockheed TV-1 
Lockheed P2V-3 
Lockheed P2V|4 
Grumman F8F-2 
Grumman F9F-2B 
McDonnell F2H-2 
Douglas AD-4 


All data in these figures were obtained during training operations unless 
otherwise stated. In order to prepare figure 24 the total number of 
load -factor peaks exceeding a given load factor N were plotted against 
the parameter ^An/An ^ 2 for each airplane. A straight line was then 
fitted to the data above a value of ^sa/An ^) 2 of about 0 . 2 , and the 
straight line extrapolated to the value N = Nq ' at zero -load -factor 
increment. The ratio n/Nq' of the number of load -fact or peaks exceeding 

a given load factor to the extrapolated value of exceeding a value of 
An/An^ =0, is plotted for each airplane in figure 24. The solid sym- 
bols in figure 24 represent less than 10 points and are considered less 
representative of the peak distributions than the other points. Also 
shown in figure 24 are the number of hours of flight represented by the 
data and the value of NoJ/T> the extrapolated value of the number of 
peaks to exceed 1 g per hour. 
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For large values of load factor the peak distribution is given in 
the form of equation ( 13 ): 



(14) 


where 


Nr 


_ t 

- r 0 T T 


and 


(15) 


w o/ T = f QY 0-6) 

When plotted on semilog paper as in figure 24, equation (l4) is a 
straight line with the slope — !■_. The values of the root mean square cr 

2a 2 - 

for the data shown vary between 0.25 and. 0.53- The average value of a 
in figure 24 is 0.284 and it appears that the straight line fits most of 
the data fairly well above la^An/An l) 2 = 0 . 0806 ^ at least up to the 

limit load factor. Of course, the value of a is to some extent a func- 
tion of the mission of the airplane and it is possible that, for a spe- 
cific mission, considerably different values than those shown in figure 24 
could be obtained. For example, in figure, 24, the training data for jbhe 
F-84E and G airplanes were obtained in Korean operations in which the 
airplane was used mostly in air-to-ground missions. The value of the 
root mean square for these data appears to be considerably lower than 
that for the other airplanes. (This effect would also be apparent if the 
true service limit load factor was lower than that assumed (7-33) in fig- 
use 24.) Most of the other training data shown in figure 24 for the 
other airplanes were obtained in training which included many types of 
missions . 


It may be noted from equation (l 6 ) that, if it can be assumed that 
fQ is about 0.08 and does not change, the percent of maneuvering flight 
time can be determined from the values of T given in figure 24. 

Tm = nq'/t = ■ V/t (17) 

T f 0 ( 0 . 08 ) ( 36 OO) 

From the values of No/f in figure 24 and from equation ( 17 ) it was 
determined that approximately 2 to 8 percent of the total flight time 
was spent in maneuvering in the Air Force training operations, about 
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2 to 6 percent was spent in maneuvering in combat operations in Korea, 
and between 7 and 24 percent in the MCA tests in operational training. 
The value of 52-percent maneuvering -flight time obtained for the five 
flights of figures 20 and 21 are applicable to Just those particular 
flights which were more active than the average. This percentage is 
obviously a function of the airplane mission and would be higher for 
gunnery or dive bombing training than for other operations. 

In addition to the data available from Air Force operations there 
are data obtained with V-G recorders from U. S. Navy operations. (See 
ref. 11.) These data are presented in figure 25 as plots of the time 
to exceed 1 g in maneuvers to ' to the time to exceed a given load 
factor t. Since, in V-G records a true frequency count cannot be 
obtained at low load factors, the data were used only where the count 
appeared to be approaching the true frequency. The times to exceed a 
given load factor were plotted against ^Zhi/dn^)^ and a straight line 
was fitted to the data and extrapolated to zero -load -factor increment 
and then the ratio of this extrapolated zero time tQ to the time to 
exceed a given load factor was plotted in figure 25 . It may be seen 
that 



( 18 ) 


The straight line shown in figure 25 is that obtained from figure 2k 
with a value of cr = 0.284. The black symbols represent points which 
consist of less than 10 peaks. The line appears to fit the data fairly 
well for most of the airplanes shown. 

In figure 26 , data obtained from references 10 and 11 are shown 
for several large airplanes. Also shown is the line obtained from the 
data for the fighter -type airplanes of figures 24 and 25 . It may be 
seen that the data for the AJ-1 and P2V-4 airplanes are close to the 
line; however, the data from the B-45 and P2V-5 airplanes appear to have 
a higher root -mean-square value than that of the line shown. It is evi- 
dent that the line shown would probably not be a good estimation of the 
loads to be encountered in airplanes which are primarily nonmaneuvering 
airplanes such as large bombers, but would be of more use in more maneu- 
verable airplanes. In addition, it is not believed that the results for 
other low -limit load-factor airplanes would necessarily indicate more 
severe use (as indicated for the B-45 and P 2 V -3 airplanes) than the line 
shown. It is believed that for any large bomber the type of mission 
would have a large effect on the slope of the line. 


Other methods of treating load- factor peak distributions .- From the 
spectrum analysis it was indicated that the major peak-load distribution 
could be given by equation ( 13 ) above lcr and this is also indicated 
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from figures 2k and 25 . It can be seen that below la, however, the 
peak-load distribution is not well represented by equation ( 13 ) . One 
way of accounting for the peaks at low load factors is to assume that 
the major peak distribution is made up of two peak distributions having 
different values of a. If such is the case, the probability of exceeding 
a load-factor peak greater than 1 g can be given as 


P 1 « 


N 


*0% 



i /An/^n L \ 2 

2 l al ) + 


1 /2xl/2*l L \ 2 

%~ 2 \ - a 2 J 

Tw 


where 


T-^ maneuvering time associated with a-j_ 

T 2 maneuvering time associated with c 2 


(19) 


For the data of figure 2k, 


a L = 0 . 28 ^ 


a 2 = 0 . 100 


T-^T = T 2 ^T = 0.50 

In figure 27 the data of figure 2k are shown with the probabilities 
given by equation ( 19 ) . The probability given by equation ( 19 ) is 
almost identical to the "Standard Probability Curve for Maneuvers" given 
in reference 1. It may be seen that the data are represented fairly 
well by- equation ( 19 ) . 


In figures 22 and 23 It might be noted that, although the major 
peaks were approximated by equation ( 13 ), the curve from the equation 
which represented the total peak count appears to be similar in shape 
to the major peak count. Therefore the major peak count could be approxi- 
mated by an equation similar to equation ( 9 ). 


The number of load -factor peaks exceeding a given load factor 
divided by the number of load -factor peaks exceeding I g from equa- 
tion ( 9 ) is 


J_ = _L _ P = 
% p o ' 1 ‘ 


1 + fn/f. 


o/ 1 


: ^M\ + f o e ' 2 a 

' N l K i* / F p 




1 - p. 


N 


) 

( 20 ) 


hi. 
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n/Nq is plotted against ^ — it is approximately linear 
2o and has a zero intercept of about ^Pq (at least for the 
of fo/^p = 0*l6o)* Therefore, if the total number of 


When 
above 

present case 

major peaks above a given load factor are plotted against Oj-J* 

and a straight line fitted to the data, the slope of the line will 
approximately determine the root -mean-square value a. Taking into 
account that the straight-line extrapolatipn of equation (20) will have 
a zero intercept of ^Pq or INq the experimental data such as those 

in figure 2k can be compared with equation (20) if it is plotted as 
N/2 Nq . In this form N/Nq = P/Pq = P]_ which is the probability of 

exceeding a given load factor for load factors greater than 1 g. In 
figure 27 the data are compared with equation (20) in this manner and 
it is in fairly good agreement. 


Thus, any of the three equations (lj), (19)^ and (20) can be used 
to approximate the major peak-load probability curve. Probably equa- 
tion ( 13 ) is the simplest if the shape of the curve below la is unim- 
portant. If the low load factors are important, for instance in some 
fatigue applications, equation ( 19 ) or (20) would be more desirable. 

In addition, probably none of these equations are valid much above the 
limit load factor since inadvertent maneuvers may affect the curves 
mostly at high load factors. (See refs. 1 and 12.) 


CONCLUSIONS 


The results of the power spectral analysis of the maneuvering load 
data Indicate that spectral methods appear to have some promise in the 
statistical treatment of maneuver loads data. It is admittedly an 
almost empirical use of these methods; however, it is an attempt to 
describe more completely the complex process of maneuvering loads on 
airplanes other than by simple peak counts . In addition, it should be 
pointed out that the flights for which the spectra were obtained were 
training flights of a general nature in which several different missions 
were simulated and that the results for an airplane flying one specific 
mission could be much different. From the present analysis the following 
conclusions have been reached. 

1. The maneuver- load-factor time histories of the five operational 
flights examined could be described by a truncated normal distribution. 

2. The power spectral densities obtained were relatively constant 
at frequencies below 0.03 cycle per second and varied, inversely with 
approximately the cube of the frequency at the higher frequencies up to 
1 cycle per second. Above a frequency of 0.005 cycle per second the 
spectrum could be approximated by a single spectrum equation. 
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3 . The frequency content above 0.2 cycle per second (5-second 
period) was very low. 

4. The load -factor peak distributions were estimated fairly well 
from the spectral analysis. 

5 . Peak -load data from about 24,000 hours of service operations 
of fighter -type airplanes appeared to agree reasonably well below the 
limit load factor with the type of equations obtained from spectrum 
analysis. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 2, 1956. 
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TABLE II.- FREQUENCY DISTRIBUTION OF TIME SPENT IN 
VARIOUS NORMAL-LOAD-FACTOR INTERVALS 


n 

Time, sec 

F-84G 

F-86A 

Flights 1 and. 2 

Flights 1, 2, and 3 

-0.2 to 

-0.001 

1 

4.5 

0 to 

.499 

63.5 

47 

.5 to 

• 599 

51.5 

29 

.6 to 

.699 

83 

70.5 

.7 to 

• 799 

124.5 

127 

.8 to 

.899 

200.5 

238.5 

.9 to 

.999 

477 

456 

1.0 to 

1.099 

1136.5 

668 

1.1 to 

1.199 

753-5 

466 

1.2 to 

1.299 

529 

332 

1.3 to 

1.399 

420.5 

274.5 

1.4 to 

1.499 

311 

239 

1.5 to- 

1.999 

987.5 

758.5 

2.0 to 

2.499 

585 

445.5 

2.5 to 

2.999 

396 

248 

3.0 to 

3-499 

275 

157.5 

3.5 to 

3-999 

117 

110 

4- . 0 to 

4.499 

90.5 

42 

4.5 to 

4.999 

36.5 

13 

5.0 to 

5-499 

18 

5-5 

5.5 to 

5-999 

10 

• 5 

6.0 to 

6.499 

4.5 


6.5 to 

6.999 

5 


7.0 to 

7.^99 

1-5 
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Figure 





Indicated airspeed, V if knots 


2.- Distribution of altitude and airspeed for flights 1 and 2 
of the F-84G airplane. 
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Pressure altitude, H p , thousands of feet 



Indicated airspeed, Vj_, knots 

- Distribution of altitude and airspeed for flights 1 , 2 , and 3 
of the F- 86 A airplane. 


Figure 3 . 
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Frequency, op« 


figure ii.- fewer spectral densities of normal load factor for three 
combined operational flights of the F-86A airplane. 
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0 *2 #4 «6 #8 1«0 

Frequency, cps 


Figure 13.- Comparison of power spectral densities of normal load factor 
on a linear scale for two operational flights of the F-84G airplane 
and three operational flights of the F-86A airpla ne . 



Frequency, ops 


Figure 14 . - Comparison of the combined power spectral densities for two 
operational flights of the F-84G airplane and three operational 
flights of the F-86A airplane with the analytical spectrum 


> ^an/ 1 + (f/r B ) a - 









Figure 15.- Illustration of truncated no rmal distribution used in 
describing threshold counts of normal load factor for F-86A air- 
plane, flights 1 , 2 , and 3- 
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(a) Flight 1. 



0 1 2 3 4 5 6 


(Id) Flight 2. 



Incremental load factor. An 


(c) Flights 1 and 2. 

Figure 16 .- Probability that a given threshold value of load factor vill 

be exceeded. F-84G airplane. 
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(a) Plight 1 . 


(h) Flight 2 . 



012345 



Incremental load factor, to 


(c) Flight 3. (d) Flights 1 , 2 , and 3. 

Figure 17.- Probability that a given threshold value of load factor will 

be exceeded. F-86A airplane. 









tti j . — -i n • < « -< • . _ _ 

rxguxe xo.- jjxsuriDuiion or normal load factor threshold counts and 
major peaks for the F-84G airplane, flights 1 and 2. 
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Figure 20.- Amount of flight time spent above given load factors for 

flights 1 and 2, F-84G airplane. 
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Figure 21.- Amount of flight time spent above given load factors for 
flights 1 , 2 , and 3 , F-86A airplane . 
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_ F-S4G- Plights 1 and 2 

O Total peak dist. | 

□ Major peak\d±st. (Method 1) 
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Normal load factor, n 


Figure 22.- Comparison of load factor peak distributions with those 
obtained from power spectral densities. F-84G airplane. 


ft , t . * 









1 2 3 4 5 6 7 


Normal load factor, n 

Figure 23 •- Comparison of load factor peak distributions with, those 
obtained from power spectral densities. F-86A airplane. 
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Figure 24 .- Percentage of load factor peaks exceeding a given value 
° f (An/ An. obtained in squadron operations for fighter-type 
airplanes. 
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O F9F-2B 

1370 

7.1 

AD-4 (Combat 
operations) 

2936 

7.i 

£2 AD-4 7 

1552 

7.1 

O F8F-2 

1092 

6.; 

OF2H-2 

5444 

6.‘ 

Q TV-1 

5660 

7.; 


•Solid symbols represent 
. peak counts of less than 10 
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Figure 25 . - Percentage of load factor peaks exceeding a given value 
° f obtained from V-g records obtained during squadron 

operations for fighter-type airplanes. 
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